Three-dimensional turbulent relative dispersion by GOY shell model 
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We study pair dispersion in a three-dimensional incompressible high Reynolds number turbulent 
flow generated by Fourier transforming the dynamics of the GOY shell model into real space. We 
show that GOY shell model can successfully reproduce both the Batchelor and the Richardson- 
Obukhov regimes of turbulent relative dispersion. We also study how the cross-over time scales 
with the initial separations of a particle pair and compare it to the prediction by Batchelor. 
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I. INTRODUCTION 



Passive particles trajectories in a fluid flow provide an 
easier way of understanding the nature of the flow than 
a full-fledged study of the Eulerian field variables of the 
flow. Recent development in Lagrangian measurement 
techniques has made it possible to study relative motion 
of pairs of fluid elements in a turbulent flow [l|, Q . The 
understanding of the pair dispersion in turbulence is of 
importance for comprehending many natural processes 
such as transport and mixing in natural and engineering 
flows, formation of warm clouds etc. 
Suppose, L is the integral scale and 77 is Kolmogorov 
length scale in a turbulent flow, then one can classify 
the process of dispersion in the flow into three distinct 
regimes based on the separation of the particles relative 
to the turbulent scales: (i) dissipation range correspond- 
ing to r{t) <^ rj {r{t) is the separation of a pair of fluid 
elements or particles at time t), (ii) inertial range corre- 
sponding to ?7 <C r{t) <g; L, and (iii) diffusion range cor- 
responding to r{t) ^ i. From the work of Richardson [3] 
(and subsequently Obukhov [4]), in the inertial range of 
fully-developed three-dimensional turbulence, there ex- 
its the well-known Richardson-Obukhov law (R-0 law): 
(r(i)^) = get^. Here, e is the energy dissipation rate per 
unit mass, and g is termed as the Richardson constant. 
Using similarity arguments, Batchelor Q extended the 
study to both short-time and intermediate-time, and ob- 
tained 
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where C2 is the Kolmogorov constant for the longitudinal 
second-order velocity structure function and 
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the classical correlation time of an eddy of size tq, is the 
time for which a particle-pair "remembers" their initial 
separation, to is also called Batchelor time. In this pa- 
per, we investigated numerically this "cross-over time" 
between the two scalings, ^ t^ and ~ t^ , and study how 
it is related to the initial separation tq. In what follows, 
we denote the first relation in equation ^ as Batchelor's 
scaling law; the second one is obviously R-0 law. Even 
though we shall be focused on three dimensional fully- 
developed turbulence, it is worth keeping in mind that 
R-0 is omnipresent [6|: there are reports of observing 
it even in two-dimensional, quasi-two-dimensional and 
rotating turbulent flows Q- Batchelor's law is no less 
robust — it has been detected even in turbulent flows 
that are not perfectly homogeneous and isotropic Q. 
In the past decade, theoretical researchers have under- 
stood — in signiflcant details — the process of advection 
of particles by developed turbulent flows in the frame- 
work of Kraichnan model[i-[il|. Therein, calculations of 
the anomalous intermittent exponents for structure func- 
tions in the inertial range show that pair structure func- 
tion, which determines the Richardson scaling exponent, 
is not intermittent. 

Until now, laboratory experiments have neither con- 
flrmed nor refuted the existence of R-0 scaling. This, 
of course, is due to the fact that these experiments are 
unable to attain very high Reynolds number due to prac- 
tical constraints. Though environmental experiments can 
achieve the highest Reynolds numbers, the problem with 
such experiments are many — inhomogeneity, anisotropy, 
rather uncontrollable experimental conditions, vastness 
of the domain to be studied, etc.. Direct numerical sim- 
ulation of Navier-Stokes equation (DNS), are much more 
successful in this respect since the computational power 
is on a rise, but cannot yet compete with the value of 
Reynolds numbers obtained in shell model simulations. 
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Also, there are issues with the commonly adopted peri- 
odic boundary conditions in DNS as they may interfere 
with the homogeneity of the system. We believe that 
shell models are well suited to test the R-0 law as not 
only can one achieve very high Reynolds number but also 
there's no sweeping effect disturbing them. It is worth 
mentioning that sweeping effect 'spresence makes it hard 
to observe Kolmogorov scaling [15| with which R-0 law 
is so intimately related being the only choice for pair dis- 
persion consistent with Kolmogorov scaling theory. A 
technique based on a real space transformations of shell 
models have recently been applied to study force fields 
study turbulent relative dispersion Kd\ and we adopt a 
similar method here. 



II. THE MODEL 

To study how a pair of particles diffuse in a 3D tur- 
bulent flow, we consider the kinematics of pair particles 
advected by the homogeneous turbulent flow obtained 
by a real-space transformation of the GOY shell model 
[l6l - [l8j . This model proposed originally by Gledzer, Ya- 
mada and Ohkitani It/, 2u\ provides a description of the 
turbulent motion embodied in the Navier-Stokes equa- 
tions. The GOY model is formulated on a A^-discrete set 
of wavenumbers, fc„ = 2", with the associated Fourier 
modes m„ evolving according to 
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Cn • En = 0, Vn. In our numerical computations we con- 
sider the vectors c„ and e„ quenched in time but averaged 
over many different realizations of these. 

As an example of the motion in this field. Fig. [T] shows 
the trajectories of two passively advected particles. As 
the relative distance diverges in time, the two particles 
experience different force fields, which in turn typically 
increase the difference in the relative velocities of the 
two particles. The figure shows the individual particles 
as they are advected, first together and later diverging 
away from each other when they are encased in different 
eddies. 

Ar-tt)=7(t)--r2(t) 
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FIG. 1: {Colour Online) Two particles being advected in a 
(3) random force field, generated by the GOY shell model. 



for n = I,-- - ,A^. The coefficients of the non-linear 
terms are constrained by two conservation laws, namely 
the total energy, E = ^„ IwnP, and the helicity (for 
3D), H = E„(-l)"fcn|W"|, or the enstrophy (for 2D} 



2D), 



Z = /c^Ju„p, in the inviscid limit, i.e. f = u = 
Therefore, they may be expressed in terms of a free pa- 
rameter only S G [0,2], a„ = 1, ^n+i = — ^, c„+2 = 
— (1 — (5). As observed by Kadanoff [22], one obtains the 
canonical value S ~ 1/2, when the helicity is conserved. 
The set ^ of A-coupled ordinary differential equa- 
tions can be numerically integrated by standard tech- 
niques [2^. We have used standard parameters in this 
paper N = 27,i^^ 10-^ ko ^ 0.05, / 5 x 10-'''(1 + i). 

The GOY model is defined in /c-space but we study 
particle dispersion in direct space obtained by an inverse 
Fourier transform [l3] of the form 
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Here the wavevectors are fc„ — fc„e„ where e„ is a unit 
vector in a random direction, for each shell n and c„ are 
unit vectors in random directions. We ensure that the 
velocity field is incompressible, V-w = 0, by constraining 



III. RESULTS 

As a main result from our simulations, we find that 
both Batchelor's and R-0 scalings are clearly observed, 
as shown in Fig. ©. Please note that we have plotted 
{\v{t) — ro]^) and not (lr(i)^ - rgj) on the y-axis. This is 
to make contact with the experimental results [1] which 
demonstrate that the correlation between the initial 
separation and the relative velocity of the pair may not 
be neglected when the flow is not perfectly homoge- 
neous. However, the GOY model — by construction — is 
a shell model of a perfectly homogeneous and isotropic 
turbulence. So, within the paradigm of shell models, 
probably it doesn't make much of a difference if one 
chooses to work with {\v{tY — Tq]). 

With no knowledge of the exact functional dependence 
of (]r(f) — ro]^) on various parameters of the flow, it is 
quite a tough task to exactly locate the cross over time, 
to, between the Batchelor regime and the R-0 regime. 
However, we can use the following rather crude but ef- 
fective method to see whether equation ^ is validated 
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ing explanation. In the case of our numerics, the Kol- 
mogorov scale, 77, is of the order of 10~^. The Reynolds 
number. Re ^ l/i^, is of the order of 10^. One may thus 
note that the Taylor microscale, A ^ Re^^^r] is more than 
100 times greater than the Kolmogorov dissipation scale. 
This leads to the fact that in our model A ~ 10~^, im- 
plying that the data points in the flat region of Fig. ^ 
do not quite correspond to the initial separations which 
are deep inside the inertial regime; and hence, it is not 
surprising that they do not follow the Batchelor predic- 
tion. Nevertheless, we feel it is quite remarkable that 
the GOY shell-model of turbulence can reproduce this 
law very well for larger separations. Similar result is ob- 
tained if one repeats the aforementioned algorithm for 
the red curve, (|r(t) — rop) /t^, in Fig. (0) focusing on 
the Batchelor's scaling regime. 

Another way to find the relation between the tq and to 
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FIG. 2: {Colour online) The initial separation ro is 10~^ and 
the data are averaged over 10000 different simulations. One 
can distinctly note the presence of both the R-0 law and the 
Batchelor's scaling law. The cross-over time to is somewhere 
in the third decade on the x-axis. The deviation of the blue 
curve (marked (|r(t) - ro|^) /t^) from R-O scaling at higher comes from the fact that R-0 scaling is (r^) = get^ while 
times is due to the fact that the pair-separations are no longer Batchelor considered {\r{t)'^ - r§ |) . So, if we plot (r^) 
in the inertial regime. All the units are in GOY-shell-model ^|^(^)2 _ j.2|^ ^^^g^g ^j^^^ ^j^^y ^-^^ become indis- 

tinguishable when R-0 scaling sets in and this will hap- 
pen around to for the corresponding tq. Thus, one can 

_ . — . — . — . . — . — . — , , , I conveniently devise an algorithm to find to in this case. 

; We may mention that also from applying this indirect 

; method, Eq. 1^ stands validated. 

In the passing, it may be mentioned that within the 
jn. _ framework adopted in this paper, one can also study exit- 

, .--*■''','** time statistics [23. [25I of the pair dispersion phenomenon 

, , , ^ •*..--■''' * to good effect. The corresponding results will be reported 

elsewhere. 
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IV. CONCLUSION 



FIG. 3: {Colour online) The cross-over times to plotted the 
initial pair-separations ro. The dashed line has the theoretical 
slope predicted by Batchelor. All the units are in GOY-shell- 
model units. 



by our simulations. To investigate where the crossover 
happens, we estimated the slope along the blue curve, 
(|r(t) - Top) /t^ in Fig. Obviously, when the slope 
is close to zero, we are in R-O scaling regime. As the 
crossover occurs before the slope is zero, we decide on the 
following algorithm: determine the time corresponding to 
the point where the slope is equal to —0.5 27| and simply 
denote that to; repeat the procedure for the curves with 
various different tq. We have used the data generated 
this way to obtain Fig. where we plot the cross-over 
times to versus the initial pair-separations tq. One can 
note the reasonable agreement of the data with the the- 
oretical Batchelor prediction (the dashed line) when the 
initial separation is not too small. What has been shown 
is that our data follow the expected trend; of course, 
data points are scattered around the theoretical slope. 
The deviation at smaller initial separations — between 
ro 



10--^ to ro = 5 X 10" 



probably has the follow- 



We have studied pair dispersion in a turbulent flow, 
applying the GOY shell model (Eq. ([3])) and Fourier 
transforming the complex shell-velocities back into real 
space (Eq. Qj). This procedure results in strongly turbu- 
lent velocity held where the dispersion of pair of particles 
can be easily studied by advecting the passive particles 
in the velocity field. In particular, we have investigated 
how the dispersion depends on the initial separations. 
We may mention that this is probably for the first time 
that a shell-model of turbulence has shown the simultane- 
ous existence of both the R-0 and the Batchelor regimes 
in turbulent pair dispersion. Whereas devising a better 
method of finding to would be worth investigating, the 
temporal extent over which these laws are valid can eas- 
ily be increased if one uses a larger number of shells. 
In closing, we hope that our results may spark interest in 
implementing GOY model to study pair dispersion in the 
enstrophy cascade dominated regime in two-dimensional 
flows. Also, it might be fruitful to use GOY model modi- 
fled appropriately 126^] for investigating turbulent relative 
dispersion in more realistic case of rotating turbulence 
which is of profound geophysical and astrophysical inter- 
est. 



4 



Acknowledgments ical results. 

We wish to thank Prof. R. Benzi for pointing out the 
importance of Taylor microscale in explaining our numer- 



A. La Porta, G.A. Voth, A.M. Crawford, J. Alexander [16] 
and E. Bodenschatz, Nature 409, 1017 (2001). 

B. Liithi and J. Berg and S. Ott and J. Mann, Physics 

of Fluids 19 045110 (2007). [17 
L.F. Richardson, Proc. R. Soc. London A, 110, 709 [18 
(1926). 

A.M. Obukhov, Izv. Akad. Nauk SSSR, Ser. Geogr. Ge- 
ofiz., 5, 453 (1941). [19 
G. K. Batchelor, Q. J. R. Meteo. Soc, 76, 133 (1950). [20 
J.P.L.C. Salazar and L.R. Collins, Annu. Rev. Fluid 
Mech. 41, 405 (2009). [21 
F. NicoUcau, G. Yu and J.C. Vassilicos, Fluid Dyn. Res. [22 
40, 68 (2008). 

N. Ouellette, H. Xu, M. Bourgoin and E. Bodenschatz, [23 
New Journal of Physics 8, 109 (2006). 

K. Gawcdzki and A. Kupiainen, Phys. Rev. Lett. 75, [24 

3834 ( 995). 

D. Bernard, K. Gawedzki and A. Kupiainen, Phys. Rev. [25 
E 54, 2564 (1996). 

L. Ts. Adzhemyan and N. V. Antonov, Phys. Rev. E 58, [26 
7381 (1998). 

N. V. Antonov and J. Honkonen, Phys. Rev. E 63, 036302 
(2001). [27 
L. Ts. Adzhemyan, N. V. Antonov, M. Hnatich and S. 
V. Novikov, Phys. Rev. E 63, 016309 (2000). 
N. V. Antonov, J. Phys. A: Math. Gen. 39, 1 (2006). 
U. Prisch, Turbulence: The Legacy of A. N. Kolmogorov, 
Cambridge University Press, Cambridge, (1995). 



M.H. Jensen, K. Sneppen and L. Angelutha, "Kol- 
mogorov Scaling from Random Force Fields" , Europhys. 
Lett, 84, 10011 (2008). 

M.H. Jensen, Phys.Rcv.Lctt. 83, 76 (1999). 

T. Bohr, M.H. Jensen, G. Paladin and A. Vulpiani, 

Dynamical Systems Approach to Turbulence, Cambridge 

University Press, Cambridge, (1998). 

E. B. Gledzer, Sov. Phys. Dokl. 18, 216 (1973). 

M. Yamada and K. Ohkitani, J. Phys. Soc. Japan 56, 

4210 (1987); Prog. Theor. Phys. 79, 1265 (1988). 

L. Bifcralc, and R.M. Kerr, Phys. Rev. E 52, 6 (1995). 

L. Kadanoff, D. Lohsc, J. Wang, and R. Benzi, Phys. 

Fluids 7, 617 (1995). 

D. Pisarenko, L. Biferale, D. Courvoisier, U. Frisch, and 
M. Vergassola, Phys. Fluids A5, 10 (1993). 
G. Boffetta and I. M. Sokolov, Phys. Rev. Lett. 88, 
094501 (2002). 

G. Boffetta, A. Cclani, A. Crisanti and A. Vulpiani, Phys. 
Rev. E 60, 6734 (1999). 

S. Chakraborty, Lev Davidovich Landau and His Impact 
on Contemporary Theoretical Physics, Horizons in World 
Physics 264, Chapter- 14, Nova PubHshers (2008). 
It may be mentioned that while experimenting with this 
approximate method of locating crossover time, we noted 
that results look more or less same for other values for 
"cut-off" slope around —0.5; thus we report herein results 
only for the cut-off slope= —0.5. 



